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Figure 1: Un-renormalized zero temperature potentials for different values of /3 and renormalized overall 
T = 0-potential (lower black points) in units of ro. The solid line shows the string potential as discussed in 
the text. 

1. Introduction 

In-medium properties of heavy quarks and heavy quark bound states are of fundamental inter- 
est for the understanding of strongly interacting matter at high temperatures and densities probed 
in current and future heavy ion collision experiments at RHIC, FAIR and LHC. The temperature 
dependence of heavy quark free energies, potential energies as well as screening masses and radii 
are important ingredients for studies on strongly coupled quark gluon plasma, possible existence 
of heavy quark bound states in the QGP as well as on transport properties of heavy quarks in the 
high temperature phase. 

We present preliminary results on heavy quark free energies §\§ that are based on an analysis of 
gauge field configurations generated by the RBC-Bielefeld collaboration in (2+l)-flavor QCD for 
the calculation of the QCD equation of state [J20. The pion mass is about 220 MeV and the strange 
quark mass is adjusted to its physical value. The calculations were performed with improved stag- 
gered fermions on lattices with temporal extent N T = 4 and 6 accompanied by high statistics zero 
temperature calculations to set the scale and to extract the zero temperature potential and corre- 
sponding renormalization constants. 

We perform a renormalization of the finite temperature heavy quark free energies and of the 
Polyakov loop using the renormalization constants obtained at zero temperature. We will dis- 
cuss the relation between heavy quark free energies, entropy and internal energy contributions 
and analyze their critical behavior in the transition region. The temperature dependence of the 
interaction between static quark anti-quark pairs will be analyzed in terms of, in general distance 
and temperature dependent, running couplings and screening masses. A preparatory study of the 
density dependence of screening masses in 2-flavor QCD with large quark masses indicates that 
non-perturbative effects in the behavior of screening is dominated by the gluonic sector. 

2. Zero temperature potential 

In fig. [I] we show the zero temperature potential, Vq q (r), in units of the distance scale ro for 
various values of the coupling j3 corresponding to different values of the lattice cut-off varying 
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Figure 2: Renormalization constants, Z len (/3 ), obtained from the matching of the zero temperature potentials 
(left). Dimensionless combinations, ro/n and ro^fa. The lower axis denotes the scale (ro/a) -1 and the 
upper axis the lattice cut-off in physical units. The lines are explained in the text (right). 



from a ~ 0.3 fm down to a ~ 0.05 fm. The scale r^/a is define by the slope of the potential, 

S^iM) =1.65, (2.1) 

and can be used to convert to physical scales with a value ro = 0.469(7) fm [[3p. Additional param- 
eters obtained from the slope of Vq q (r) are the string tension, oa 2 , and r\/a defined by replacing 



1.65 by 1.0 in ( |2.1| ). The dimensionless combinations of these parameters displayed in fig. |2j show 
only small cut-off effect, e.g. roy^ stays constant in the entire range of couplings in which the 
lattice spacing changes by a factor of 6. Using a quadratic fit Ansatz for a < 0.1 5fm we obtain 
rQy/a = 1.1034(40) and ro/n = 1.4636(60). For both ratios we observe that 0(a 2 ) corrections 
are small. 

The potentials as calculated on the lattice are ultraviolet divergent and need to be renormalized. We 
have matched all potentials to a common value at large distances, r/ro = 1.5, taken to be identical 
to the large distance string potential which in units of ro is given by 

roV string (r/r ) = -rr^-j- + (or 2 )-, (2.2) 
12r/r r 

where we have used the value roy^ quoted above. The result of this renormalization is shown 
in the lower part of fig. |T} The good matching of all potential data over the whole distance range 
again shows that cut-off effects are small in this observable. 

The matching procedure provides renormalization constants for the static quark potential which we 
will use for the renormalization of Polyakov loops and heavy quark free energies at finite tempera- 
tures. They are defined as Z ren (/3) = exp((c(j8)a)/2), where c(j8) denotes the constant shift for the 
corresponding zero temperature potential. Already note here that no additional divergencies arise 
when going from zero to finite temperatures. The solid line in fig. [j] shows the string potential. 
Note that the deviations of the potential V qq (r) from V str m g (r) at small distances shows the effect 
of the QCD running coupling while the Coulombic term in the string potential stems from large 
distance corrections in the string model. Furthermore we do not observe string breaking here which 
is due the operator used in our analysis and the limited distance range analyzed here. 
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Figure 3: Renormalized heavy quark free energies, F\ (r, T), for different values of the temperature (left). 
The solid line shows the zero temperature potential V\(r). Asymptotic values for the heavy quark free 
energies (right), F\(r = °°,T), for (2+1) flavors compared to earlier results for pure gauge theory and 2- 
flavor[|||. 

3. Heavy quark free energies 

The heavy quark free energy of a static quark-antiquark pair in a color singlet state separated 
by distance r is defined by 

= -log (Tr (L ren (0)L ren (r))) , (3.1) 
where L ren is the renormalized Polyakov loop, 

N T -l 

L ren (f) = (ZrenO?))^ L bare = ]J Z ma {$)U {1Mfl . (3.2) 

-vo=0 

The operator used in ( fO| ) in general is gauge dependent and we fix to Coulomb gauge accord- 
ing to arguments based on [[7], g]. As already noted we are using the renormalization constants, 
Zren(/3), obtained at zero temperature for the renormalization of the Polyakov loop. In this way 
also the heavy quark free energies are properly renormalized, which is evident from fig. [| (left) 
where F\ (r, T) for different temperatures is shown. The solid line is the zero temperature potential 
V qq (r). 

The free energies become temperature independent at small separations and coincide with Vq q (r) 
showing the correct renormalization. At larger distances temperature effects set in, below T c due 
to string breaking and above T c due to screening of the static sources in the thermal deconfined 
medium. The onset of this temperature effects is shifted towards smaller distances with increasing 
temperature. The asymptotic values, F\(r = °°,T), are show in fig. ^] (right) compared to previous 
results from quenched and 2-flavor QCD with larger quark masses ||, ||]. The qualitative be- 
havior is comparable in all those theories. Note that F^ is infinite for the quenched theory below 
the critical temperature. With dynamical quarks below the critical temperature F^ is already close 
to the value estimated from zero temperature (indicated by the black band), it shows a strong de- 
crease around T c and a quite linear behavior at high temperatures. This temperature dependence 
already indicates that entropy contributions play an important role especially around the Uansition 
temperature. 
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Figure 4: The renormalized Polyakov loop for different temporal lattice extents, N T = 4,6 and 8 (left). The 
vertical lines show the location of the transition temperature determined in |Q| on lattices with temporal 
extent N T = 4 (right line) and N T = 6 (left line), respectively. Comparison to earlier results (right) for pure 
gauge theory and 2-flavor [Q, ||, |(|. 



4. Renormalized Polyakov loop 

As already discussed in the previous section, the renormalization constants obtained at zero 
temperature can be used for the renormalization of the Polyakov loop at finite temperature ac- 



cording to (p.2p. Due to the cluster property this renormalization is equivalent to a definition 
of L ren using the asymptotic value of the heavy quark free energy, F^T) = F\(r = °°,T), i.e. 
^ren = ex p( — F(r = °°,T)/2T). While L ren is zero below the critical temperature in pure SU(3) 
gauge theory and has a finite gap at T c due to the first order transition in this theory, it is non-zero 
in QCD with dynamical quarks even in the low temperature phase. Although the temperature de- 
pendence is continuous here, L ren shows a pronounced rise around the transition region from small 
to large values in the high temperature phase. 

We note that the most rapid change is in good agreement with the region where the chiral con- 
densate as well as bulk thermodynamic quantities, e.g. energy and entropy densities change most 
rapidly Furthermore the cut-off dependence of L ren on lattices with temporal extent N x = 4 
and 6 is small, which is in agreement with results obtained in studies of L ren in pure SU(3) gauge 
theories [||,|9|]. The large cut-off dependence observed in a study with the 1-link stout smeared stag- 
gered action used in JTOj ] mainly seems to arise from the cut-off dependence of the zero temperature 
observable (fg) used to set the temperature scale. The fact that the renormalized Polyakov loop 
becomes larger than one at high temperatures was already observed in the quenched theory and 
is predicted by perturbation theory where the high temperature limit, L ren (T = °°) = 1, is reached 



from above [ |11| , |12Q . Note that the renormalized Polyakov loop as defined in ( p.2| ) is no longer a 
SU(3) matrix. 



5. Entropy and internal energy contributions 

The temperature dependence of the heavy quark-antiquark free energy at asymptotic large dis- 
tances, Fco(T), already indicated that entropy contributions play an important role especially in the 
vicinity of the transition region. Thus it becomes important for discussions on strongly coupled 



5 



Screening at finite temperature and density 



Olaf Kaczmarek 



5000 
4500 
4000 
3500 
3000 
2500 
2000 
1500 
1000 
500 


-500 



UJMeV] 
FJMeV] 



4000 
3500 
3000 
2500 
2000 
1500 
1000 
500 


-500 



TS,.[Me\»] 
F.,[MeV] 







Figure 5: Asymptotic large distance values for internal energy (left) and entropy (right) compared to heavy 
quark free energies. The left band indicates the zero temperature limit of the energy and the red line shows 
the vanishing T = limit for TSo. 



quark gluon plasma (SQGP) [ |13| , |14| , |15Q , the possible existence of heavy quark bound states that 
might survive above deconfinement [16, 17|, [l8|,[[9|, 20, |l]] as well as on transport properties [22] 
in the high temperature phase. 

Here we will only discuss the asymptotic large distance behavior of these quantities, i.e. Soo(T) = 
Si (r = oo 5 T) and U^(T) = U\{r = oo ; T). To separate entropy, S^, and internal energy, U^, contri- 
butions from heavy quark free energies we use standard thermodynamic relations, 



SJT) 



dF„o(r) 



and Uo.iT) 



dF^{T)/T 



(5.1) 



dT v 7 dT ' 

to calculated their asymptotic behavior. In fig. |5] we compare internal energy (left) and entropy 
(right) to the asymptotic behavior of the free energy. The band at the left axis indicates the value of 
the energy in the zero temperature limit and the line in the right figure indicates that entropy, TS^, 
vanishes at zero temperature. 

In contrast to the monotonic decrease observed for Foo, both £/„ and TSoo show qualitative different 
(critical) behavior in the transition region. While at small temperatures they rapidly approach 
values close to their zero temperature limits, both show a pronounced peak, Uoc(T ~ T c ) ~ 4.5 GeV 
and TSoo(T ~ T c ) ~ 4 GeV, which again decreases rapidly toward smaller values and become rather 
flat above 1.5 T c . Note that the temperature where both observables attain their maximum is in 
good agreement with the critical temperature obtained using different observables in [§]. Similar 
behavior was already observed in previous studies with «/=2 and 3 with larger quark masses in 



,23]. 



6. Effective running coupling constant 

As already discussed, the heavy quark-antiquark free energies become temperature indepen- 
dent at small separations. To analyze the onset of temperature effects it is more appropriate to study 
the effective distance and temperature dependent running coupling defined through 



0£eff(r,r) 



3 2 dF l (r,T) 
—r 

4 dr 



(6.1) 
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Figure 6: The effective (distance dependent) running coupling (left). The solid line shows the zero tem- 
perature running coupling calculated from the zero temperature potential. Effective temperature dependent 
coupling (right), a max {T), define by the maximum of G!eff(r, T) as explained in the text for (2+1) flavor com- 
pared to 2-flavor and pure gauge theory [[|, |], 0]. The lines show a perturbative inspired fit-Ansatz for high 
temperatures. 

In fig. H (left) results for various temperatures are compared to the zero temperature running cou- 
pling, OCT=o(r), defined in the same way (solid line). Note that the quadratic rise of (Xt=o is a 
non-perturbative effect that stems from the linear rising string tension term in the potential, while 
at small distances the logarithmic weakening of the coupling is visible and at sufficiently small 
distances it should reach the perturbative (asymptotic free) behavior. 

At finite temperatures, a e ff (r, T) follows this zero temperature behavior to relatively large distances, 
before screening sets in leading to a maximum and a decrease at larger distances. Although the on- 
set of temperature effects and the maximum is shifted towards smaller distances with increasing 
temperature, at temperatures slightly above the critical one, a e ff(r, T) still follows the quadratic 
behavior quite far, indicating that even above T c remnants of confinement forces are present. 
We have used the value of the maximum of the running coupling, a max (r), to define an effective 
temperature dependent coupling constant that can be used to indicate an effective coupling strength 
at distances where the screening of the static quark-antiquark pair sets in. The results are shown 
in fig. ^ (right) in comparison with pure SU(3) gauge theory and 2-flavor QCD results from [j|]. 
The solid lines show results of a fit with an Ansatz using the two-loop perturbative coupling at 
temperatures T > l.2T c . Note that the 2-flavor results are calculated for rather large quark masses 
(m/T = 0.4). Therefore the increase of (Vax for (2+l)-flavor is rather a quark mass effect than 
a flavor effect. The comparison of the results for temporal extent N t = 4 and 6 again shows that 
cut-off effects also in this observable are quite small. The large value of the coupling close to the 
critical temperature should not be identified with a large Coulombic coupling but clearly includes 
non-perturbative effects. 

7. Debye screening at finite temperature and density 

We follow the commonly used approach to define the non-perturbative (Debye) screening 
mass, mj)(T), and temperature dependent coupling constant, a(T), using a screened Coulomb fit 
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Figure 7: Debye screening masses in units of the temperature (left) for different number of flavors. The 
solid lines and error band show the fit of a perturbative Ansatz including a scale factor A as described in the 
text. Effective T-dependent Coulombic couplings (right) as obtained from the same fits compared to (Xmax- 

Ansatz for the large distance part of the heavy quark free energies normalized by its asymptotic 
large distance value, 

Fi (r, T) -F 1 (r = », T) = - ~ ^r^. (7.1) 

Our results for the screening masses are summarized in fig. [7J (left) as function of T /T c and are 
compared to results in pure gauge theory and 2-flavor QCD from 

Although we are not expecting perturbation theory to hold at the small temperatures analyzed here, 
the enhancement for increasing number of flavors is in qualitative agreement with leading order 
perturbation theory, i.e. 

Here we have already introduced a multiplicative constant A to allow for non-perturbative contri- 
butions (A = 1 in perturbation theory). Using a two-loop perturbative definition of g(T) a best fit 
analysis for mo(T)/T for temperatures T > l.2T c leads to A = 1.52(2), 1,42(2) and 1.66(2) for 
Np = 0,2 and (2 + 1), respectively. The results including an error band are shown by the solid lines 
in fig. ^ (left). 

The large value of A indicates that non-perturbative contributions are important in the temperature 
range analyzed here. Note that a study in pure gauge theory up to temperatures as high as 24 T c led 
to an only slightly smaller value, A = 1.39(2), revealing that even at such high temperatures the 
screening mass still is far from being perturbative. 

The results for the coupling, oc(T), are shown in fig. [7] (right) in comparison with a max as discussed 
in the previous section. The large values of Ct(T), especially close to the transition region, de- 
scribes the large distance behavior of heavy quark free energies and should not be confused with 
the coupling a raax that characterizes the short distance part of F\ (r, T). The latter is almost temper- 
ature independent and can to some extent be described by the zero temperature coupling. 

So far we have only discussed results for vanishing baryon number densities. In leading order per- 
turbation theory, the dependence of the Debye mass on temperature and quark chemical potential 
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Figure 8: Zero (left) and second (right) order density coefficients of the Taylor expansion of heavy quark 
free energies in 2-flavor QCD j24|]. The solid line in the left figure shows the fit and in the right figure the 
leading order perturbative result as explained in the text. 



is given by 

Here we will discuss results for two-flavor QCD with large quark masses [Qj. The extension to 
non-zero densities in terms of a Taylor expansion in the quark chemical potential is straightforward 
and comparable to the Taylor expansion for bulk thermodynamic observables [|25|], 

m D {T^ q ) = m (T) +m 2 (T) (^)' + 0{)l% (7.4) 

where mo(T) is the zero-order (zero density) contribution as already discussed above and rri2(T) is 
the second order Taylor coefficient. Note that for a qq-system the odd order coefficients vanish. 
The results for mo (r) andra2(r) are shown in fig. [8] (left) and (right), respectively. The solid line in 
the left figure shows the fit result for mo as discussed above and the line in the right figure shows the 
second order expansion coefficient from the perturbative leading order expression (f7j|). While the 
zero order coefficient clearly shows non-perturbative behavior up to high temperatures, the good 
agreement of the second order coefficient with the perturbative expectation at temperatures above 
1.5r c indicates that non-perturbative contributions to screening in the high temperature phase of 
QCD are dominated by the gluonic sector. 

The density dependence of screening masses so far was only analyzed for 2-flavor QCD and large 
quark masses. For a more realistic description the analysis will be carried out for (2+l)-flavor QCD 
with almost realistic quark masses in the future. 



8. Conclusions 

We have analyzed heavy quark free energies, their temperature dependence and screening 
properties for (2+l)-flavor QCD with almost realistic quark masses. The renormalization of the 
free energies as well as the Polyakov loop was performed using renormalization constants obtained 
at zero temperatures. For the asymptotic large distance behavior, entropy and internal energy con- 
tributions were separated from the free energies and show critical behavior around the transition 
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region. 

From the analysis of the effective running coupling constants the onset of medium effects on a 
heavy quark-antiquark pair was analyzed and from the large distance behavior of the free energies 
we extracted screening properties at vanishing density. The extension to finite baryon densities so 
far is limited to the case of 2-flavor QCD and large quarks masses. The agreement of the second 
order expansion coefficient with perturbation theory at temperatures above \.5T C indicates that the 
main non-perturbative contributions to screening stems from the gluonic sector while fermionic 
contributions seem to be small. This behavior has to be confirmed for (2+l)-flavor QCD with 
almost realistic quark masses in the future. 

References 

[1] RBC-Bielefeld collaboration, in preparation. 

[2] M. Cheng et al. [RBC-Bielefeld collaboration], Phys. Rev. D 74, 054507 (2006) and arXiv:0710.0354. 

[3] A. Gray et al., Phys. Rev. D 72, 094507 (2005). 

[4] O. Kaczmarek et al., Prog. Theor. Phys. Suppl. 153 (2004) 287. 

[5] O. Kaczmarek and F. Zantow, Phys. Rev. D 71, 1 14510 (2005). 

[6] O. Kaczmarek and F. Zantow, hep-lat/0506019. 

[7] O. Philipsen, Phys. Lett. B 535, 138 (2002). 

[8] O. Jahn and O. Philipsen, Phys. Rev. D 70, 074504 (2004). 

[9] O. Kaczmarek, F. Karsch, P. Petreczky and F. Zantow, Phys. Lett. B 543, 41 (2002). 

[10] Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, Phys. Lett. B 643, 46 (2006). 

[11] E. Gava and R. Jengo, Phys. Lett. B 105 (1981) 285. 

[12] P. Petreczky, Eur. Phys. J. C 43, 51 (2005) 

[13] J. Liao andE. Shuryak, arXiv:0706.4465. 

[14] E. Shuryak, Prog. Part. Nucl. Phys. 53, 273 (2004). 

[15] M. Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 (2005). 

[16] S. Digal, P. Petreczky and H. Satz, Phys. Lett. B 514, 57 (2001), Phys. Rev. D 64, 094015 (2001) 

[17] A. Mocsy and P. Petreczky, Phys. Rev. D 73, 074007 (2006) and arXiv:0705.2559, arXiv:0706.2183. 

[18] D. Cabrera and R. Rapp, hep-ph/061 1 134. 

[19] Y. Park, K. I. Kim, T. Song, S. H. Lee and C. Y. Wong, arXiv:0704.3770. 

[20] W. M. Alberico, A. Beraudo, A. De Pace and A. Molinari, arXiv:0706.2846. 

[21] E. V. Shuryak and I. Zahed, Phys. Rev. D 70, 054507 (2004). 

[22] H. van Hees, M. Mannarelli, V. Greco and R. Rapp, arXiv:0709.2884. 

[23] P. Petreczky and K. Petrov, Phys. Rev. D70 (2004) 054503. 

[24] M. Doring, S. Ejiri, O. Kaczmarek, F. Karsch and E. Laermann, Eur. Phys. J. C 46, 179 (2006). 

[25] C. R. Allton et al, Phys. Rev. D 71, 054508 (2005). 



10 



